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INTRODUCTION 

Due  to  the  difference  in  the  effect  of  applied  pressure  on 
melting  points,  a  substance  which  melts  at  a  temperature  lower 
than  the  freezing  point  of  an  aqueous  solution  may  be  made  to  melt 
at  a  temperature  higher  than  the  melting  point  of  water  by  a  suf- 
ficiently high  applied  pressure.   Therefore,  a  working  medium  can 
be  selected  to  form  a  cyclic  auxiliary  system  for  removing  the 
heat  from  the  salt  water  to  produce  ice  and  concentrated  brine  and 
for  supplying  heat  to  melt  the  ice  so  formed  at  a  higher  pressure. 
At  the  end  of  a  cycle  the  auxiliary  system  is  returned  to  its 
original  state.   This  is  the  basic  underlying  theory  of  the  de- 
salination freezing  process  (hereafter,  it  will  be  called  the  In- 
version  Freezing  Process)  investigated  here. 

The  specific  objectives  of  this  research  were: 

1.  To  design  and  build  experimental  equipment  which  could 
be  used  to  produce  ice  continuously  under  atmospheric 
pressure. 

2.  To  determine  the  ice  production  rate  under  various  operat- 
ing conditions. 

3.  To  determine  the  ice  crystal  sizes  produced  from  this 
process. 


REVIEW  0?  LITERATURE 

The  oceans  represent  the  largest  water  reservoir  of  the 
earth.   Only  29%   of  the  global  surface  is  covered  with  land;  the 
remaining  7~L%   is  ocean.   Ocean  water,  however,  contains  on  the 
average  3~2  Per  cent  (by  weight)  of  dissolved  salts  mainly  NaCl,  a 
concentration  that  makes  the  water  unsuitable  for  many  household 
and  industrial  uses  and  for  the  irrigation  of  conventional  land- 
grown  crops. 

There  is  at  all  times  an  exchange  of  water  between  oceans 
and  seas,  on  the  one  hand;  and  rivers  and  lakes  on  the  other  hand, 
the  emptying  of  a  river  into  the  sea  is  an  illustration  of  the 
flow  of  fresh  water  into  salt  water.   Flow  in  the  opposite  direc- 
tion also  takes  place;  moisture  evaporates  as  a  result  of  solar 
hea     from  the  surface  of  the  sea  and  is  carried  upward  t:o  t. 
cooler  regions  of  the  atmosphere,  where  it  is  condense-  and  retu 
ed  :o  earth  again  as  rain  to  fill  the  streams,  rivers,  and  lakes 
rith  fresh  water. 

The  natural  evaporation- condensation  action  is,  by  far,  t) 
-  iportant  of  nature's  desalination  processes.   Alt        it- 
il  desalination  processes  such  as  freezing  and  biological  pro* 
.'j,  exist. 

fortunately,  fresh  water  produced  from  the  sea  ... 
Is  is  not  always  available  at  the  locations 
:  required  by  man.   In  the  past,  decisions  on 

ive,  to  a  large  extent,  been  dependent  upon  the 
.atural  supply  of  fr« -%. 


In  recent  decades,  as  the  result  of  the  large  increase  in 
world  population  and  various  military  and  political  factors, 
there  has  been  a  tremendous  drive  to  open  arid  areas  for  large- 
scale  settlement.   Moreover,  water  shortages  appear  occasionally 
in  highly  industrialized  areas  in  zones  of  moderate  climate  with 
abundant  rainfall.  Wnen   there  is  no  fresh  water  within  a  reason- 
able distance  and  ocean  water  or  brackish  ground  water  is  avai- 
lable, salt  removal  from  water  becomes  necessary.   This  process  is 
called  "salt  water  conversion"  or  "water  desalination". 

Salt  water  conversion  can  be  accomplished  in  many  ways.   Some 
methods  of  producing  fresh  water  from  salt  water  have  been  known 
in  concept  for  centuries;  other  methods  have  been  uncovered  only 
in  the  past  few  years.   Some  methods  are  now  utilized  commercially, 
some  have  progressed  to  the  pilot  plant  stage  of  development,  and 
others  are  still  under  study  in  laboratories. 

The  principal  desalination  processes  can  be  summarized  as 
follows  (1). 

A.   Distillation: 

a.  Long-tube  vertical 

b.  Flash  distillation 

c.  Vapor  compression 

d.  Mult ief feet  multistage 
5.   Membranes: 

a.  Reverse  osmosis 

b.  Slectrodialysis 

c.  Transport  depletion 
C.   Freezing: 


a.  Direct  freezing 

b.  Indirect  freezing 
D.  Humidificationi 

a.  Solar 

b.  Diffusion 
E«   Chemical i 

a.  Hydrates 

b.  Ion  exchange 

'  e  cc::72::tiomal  freezing  processes. 

Separation  of  fresh  water  from  salt  water  by  freezing  is 
based  on  the  fact  that  when  salt  water  freezes,  fresh  water  ice 
crystals  form,  and  the  salt  remains  in  solution  in  the  unfrozen 
water.   A  freezing  process  is  generally  a  two  step  operation  re- 
cjuiring  both  cooling  and  heating.   In  any  freezing  process,  th 
salt  water  is  cooled  until  ice  is  formed.   The  ice  is  separated 
from  the  brine  and  then  melted  to  produce  the  product  water  of  t 
process. 

Freezing  processes  have  been  studied  for  over  20  years  and 
are  still  in  a  period  of  rapid  development.   Generally,  they  c 
be  divided  into  two  categories,  namely,  direct  freezing  ar.d  Indi- 
rect freezing.   In  the  former,  water  acts  as  its  own  refri     z , 

,    the  evaporation  of  water  vapor  removes  th 
for   .       ition  of  ice.   In  the  latter,  a  liquid  . 
immiscible  with  water,  such  as  butane,  is 
it. 


In  both  direct  freezing  and  indirect  freezing,  ice  crystals 
are  formed  in  a  freezing  vessel  (or  crystallizer)  in  which  the 
pressure  is  maintained  at  about  0.005  at  a.   As  rapid  evaporation 
of  the  refrigerant  takes  place  in  it,  the  sea  water  cools  and 
eventually  freezes  (The  freezing  temperature  is  -2.05°C  for  pure 
sea  water  and  -3.8  C  for  doubly  concentrated  sea  water,  which  is 
the  liquid  phase  in  the  crystallizer).   In  other  words,  the  inlet 
stream(s)  splits  into  vapor  and  ice,  and  the  heat  needed  for  the 
evaporation  is  supplied  by  conversion  of  x^ater  to  ice. 

The  ice  produced  from  the  crystallizer  must  be  washed  be- 
cause it  is  in  the  form  of  a  fine  slush  still  holding  up  to  50 
weight  per  cent  brine  in  the  space  between  the  crystals.   Coun- 
tercurrent  washing  procedures  have  been  developed  in  which  the 
ice  crystals  rise  against  a  stream  of  wash  water  flowing  downward. 
This  method  is  so  efficient  that  only  a  few  percent  of  the  fresh 
water  has  to  be  used  for  this  purpose  (2). 

The  pressure  in  the  crystallizer  is  maintained  by  the  con- 
tinuous pumping  of  vapor  by  a  compressor.   For  the  purpose  of 
saving  energy  and  equipment,  the  vapor  is  pumped  into  a  vessel, 
called  a  melter,  in  which  the  pressure  is  maintained  at  4.6  mm 
mercury  in  the  direct  freezing  process  and  a  little  bit  over  atmos- 
pheric pressure  in  the  indirect  freezing .process.   The  washed  ice 
is  also  continuously  fed.   Into  the  melter  within  the  vessel,  the 
refrigerant  vapor  condenses  to  liquid;  at  the  same  time,  the  ice 
crystals  melt  and  turn  into  water.   In  the  direct  freezing  process, 
liquid  water  is  obtained  from  the  condensation  of  vapor  and  be- 
comes a  part  of  the  fresh  water  product.   In  the  indirect  freez- 


•orocess,  vaporous  refrigerant  is  condensed  and  returned  to 
the  crystallizer. 

Freezing  processes  have  the  following  characteristics. 
la   They  operate  at  low  temperatures  thus  minimizing  scale 
and  corrosion  problems.   Low  temperature  plants  may  be 
expected  to  last  longer  than  those  operating  at  hi 
temperatures. 

2.  In  the  freezing  and  melting  steps  heat  transfer  is  accom- 
plished without  passing  through  metal  surfaces.   The 
great  advantages  of  transferring  the  latent  heat  of  phase 
change  by  direct  contact  between  phases  are  that:  scale 
formation  is  unimportant  and  that  temperature  differen- 
ces required  for  heat  transfer  are  small. 

3.  Freezing  processes  can  be  designed  for  practically  any 
feed  from  dilute  brackish  water  to  sea  water. 

4.  Freezing  processes  operating  on  brackish  water  can  cc  - 
vert  economically  a  larger  fraction  to  fresh  water  pro- 
duct than  other  processes,  such  as  electrodialysis  a 
reverse  osmosis. 

5.  Freezing  is  the  only  developed  process  which  can  w 
electrical  power  directly  to  desalt  .     .ter. 
allows  plants  to  be  built  where  cheap  steam  is  not 

le.   Electrical  power  may  be  transport  o:" 

-les,  but  ste     1st  be  used  near  its  sourc  . 
Due  to  the  n      of  the  equipment  used  In 
ltable  size  of  a         process  v. 
be      50 1OOC  lions  of 


day  (That  of  a  distillation  process  appears  to  be  on  the 
order  of   100  million  gallons  per  day). 
The  conventional  freezing  processes  have  the  following  dis- 
advantages . 

1.  High  equipment  cost  due  to  the  large  volume  of  vapor 
involved. 

2.  High  energy  cost  because  of  the  low  energy  efficiency  of 
the  compressor  in  compressing  the  vapor. 

3.  Poor  control  of  crystal  formation  because  of  the  vigo- 
rous agitation  and  the  local  subcooling  due  to  the  ra- 
pid vaporization  of  liquid. 

THE  INVERSION  FREEZING  PROCESS. 

This  process,  which  was  proposed  by  the  Chengs  (3)»  appears 
to  possess  the  same  advantages  as  the  conventional  freezing  pro- 
cesses but  not  their  disadvantages.   The  process  is  based  on  a 
distinctly  different  way  of  upgrading  heat.   Due  to  the  difference 
in  the  effect  of  applied  pressure  on  melting  points,  a  substance 
which  melts  at  a  temperature  lower  than  the  freezing  point  of  an 
aqueous  solution  may  melt  at  a  temperature  higher  than  the  melt- 
ing point  of  water  at  a  sufficiently  high  applied  pressure.   There- 
fore, a  working  medium  (or  refrigerant)  can  be  selected  to  form  a 
cyclic  auxiliary  system  to  remove  the  heat  of  crystallization  in 
the  partial  freezing  operation  and  to  supply  the  heat  required  to 
t  the  ice  so  formed.   The  process  is  distinct  from  the  cor^veT.- 
.onal  freezing  processes  in  that  it  handles  only  the  condensed 
[liquid  and  solid)  phases.   This  has  a  considerable  effect  on  zhe 


: 


energy  requirments  of  the  process  and  favors  the  control  of  ice 
cry stall! zat  i  on . 

-c:-;r.   The  shape  and  direction  of  a  univariant  pressure  v  . 
temperature  curves  for  the  melting  point  of  a  pure  substance  a. 
the  eutectic  temperature  of  a  binary  mixture  are  given  by  the 
Clausius-Clapeyron  equation  (^).   The  melting  point  is  seen  to  be 
raised  (positive  slope)  as  pressure  increases  if  the  solid  is  cens- 
er than  the  liquid;  a  negative  slope  is  found  in  a  few  cases, 
such  as  water,  bismuth  and  gallium,  where  the  liquid  is  denser 
than  the  solid.   According  to  P.  W.  Bridgman  (5),  the  melting 
point  of  water  decreases  one  degree  centigrade  with  an  increase 
in  the  applied  pressure  of  about  100  atmospheres,  and  the  general 
tendency  of  most  systems  is  to  have  a  positive  slope  of  the  fus 
rve,  the  average  value  being  50-60  atm./deg.  C. 

The  Inversion  Freezing  Process  employs  a  substance  or  a  mix- 
ture of  substances  as  a  working  medium.   It  is  so  selected  that 
its  melting  point  is  lower  than  the  freezing  point  of 
solution  from  which  fresh  water  is  to  be  separated 

However,  when  the  applied  pressure  exceeds  a  cert   .       . 
Lelting  point  is  raised  to  exceed  the  melting  point  of 
he  same  applied  pressure.   Figure  1  cor. 
.nt  change  of  water  with  pressure  with  that  of  an  : 

It  also  shows  schematically  the  conditions  o  ' 

.  .re  ice-making  operation,  the  lnv     n  poir. .,  .. 
-  .elting  operation.   It  reveals  that 
of  substance  A  is  lower  th  oint  c 

If  solid  A  is  contact.-  with  liquid  B,  direct.  , 
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Phase    equilibrium   diagram  for  an   aqueous  system 
and    an  idea!    working   medium  (  schematic). 
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solid.  A  will  melt  and  liquid  3  will  be  solidified  cy   the 
exchanged  between  them.  Therefore,  a  substance  or  mixture  of  sub- 
stances which  forms  a  working  medium  oo  absorb  the  heat  of  crystal- 
lization at  a  low  pressure  by  melting  can  remove  the  heat  cz"   fusion 
of  ice  at  a  sufficiently  high  pressure  by  solidifying     If. 

Figure  2  illustrates  the  general  process  and  shows  z'r.cz   par- 
tial freezing  of  an  aqueous  solution  can  be  coupled  with  meltii 
oC   the  working  medium  by  the  following  procedure: 

Aqueous  Solution  Ice 

->-        (at  low  pressure V  +  Cone.  Aqueous  Solution 
Working  Medium  +  Working  Medium 


(solid)  (liquid) 


(1) 


melting  of  ice  is  coupled  with  the  solidification  of  the 
working  medium  through  the  following  procedure: 

Ice  Fresh  Water 

+      (at  hi"rh  pressure)         +  (2) 

Working  Medium  ~~^   Working  Medium 

(liquid)  (solid) 

Considering  the  overall  effect,  the  working  medium  oper  -     .  a 
cycle  to  reuse  heat  by  drawing  the  heat  of  crystallizatioi 

:rlng  it  to  the  melting  of  ice.   Some  wor    s  to  be  su  - 
plied  to  the  system  in  order  to  conform  to  thermod; 

2H  1        i  wry,-        j     The  work  in, 

ire  substance  or  i      ir  of  substances.  re 

siting  point,  such  i cut ec tic  point,  or  I 

point  ran   .   It  Is  desi.    i  to  sel   .am.         .;e 
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melting-point  range  approximates  closely  the  range  of  the  freezin 
point  of  the  aqueous  solution  during  the  freezing  operatic... 

An  ideal  working  medium  should  have  the  following  qualitiesi 
a  proper  melting  range,  a  very  low  solubility  in  water,  a  1 

Lue  for  the  latent  heat  of  fusion,  and  a  low  value  cf  [%m)    for 
Ltingi   It  should  also  be  non-toxic,  cheap,  and  readily  availa- 
ble.  Theoretically,  any  substance  which  has  a  proper  melting 
point  (or  melting  range)  and  has  low  solubility  in  water  can  be 
used.    it  the  problem  of  economy  will  probably  limit  the  practi- 
cal working  medium  to  a  hydrocarbon  or  mixture  of  hydrocarbons • 

Phase  diagrams  (solid-liquid  equilibrium)  for  binary  systems 
of  saturated  hycrocarbons  have  been  reported  by  R.  Salzgeber  (6) 
and  by  S.  3.  Maa  (7).   Figures  A-l,  A- 2  and  A- 3  in  Ap]  ndix  A  chow 
se  diagrams,  obtained  under  atmospheric  pressure  and  hi  . 
.2 cures,  of  n-tridecane  (pure  grade)— n-tetrade cane  (pure 
n-tridecane  (technical  grade)— n-tetrade cane  (technical  grade), 

id  n-tridecane  (technical  grade) — n-pentadecane  (technical  _rc.de) 
binary  systems,  respectively  (7).   It  can  be  seen  that  ...     ■   of 

these  systems  with  proper  composition  meets  the  required  mell 

range  both  at  atmospheric  pressure  and  high  pressures.   1 
it  of   fusion  of  high  n-paraffin  hydrocarbons  is  particu! 
.  around  50-60  cal/gr  for  even  C-number  hydrocarbo]  i 
round  kO   cal/gr  for  odd  C-number  hydrocarbons  (3).  -c- 

have  also  been  reported  to  be  insoluble  in 

curors  have  estimated  the  long  ran  - 

lydrocarbons  for  quantities  of     ral  mill i or. 
ir  co  be  within  k.5-5  t/Va   (10,11). 
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Benzene  forms  eutectic  mixtures  with  many  substances  (12), 
and  the  eutectic  temperatures  are  in  the  range  suitable  for  the 

salination  of  sea  water.   A  mixture  of  benzene  and  naphthalene 

(3  p 
is  a  desirable  system  because  it  has  a  lex  value  of  ( — - )    for 

d? 

and  is  readily  available.   There  are  objections  to  its  use 
as  a  working  medium,  however.   They  are  that  it  has  a  slight  but 
signigicant  solubility  in  water,  it  is  toxic,  and  it  has  rather 
high  volatility.   Using  this  mixture  as  a  working  medium  would  re- 
quire post  treatments  to  remove  it  both  from  the  product  fresh 
water  and  from  the  reject  salt  water. 

A  working  medium  may  be  found  among  alkylnaphthalenes ,  for 
example,  2-n-propylnaphthalene  (melting  point:  -3°  C,  density i 
0.9770  gr/ml,  very  slight  solubility  in  water).   Since  the  den- 
sity of  an  alkylnaphthalene  is  so  close  to  that  of  an  aqueous 
solution,  elaborate  equipment  may  be  required  for  their  sepa- 
r£o i  en . 

Process  description.   Figures  3  and  ^  are  the  block  diagram 
and  flow-sheet  respectively  of  the  proposed  process.   The  process 
is  divided  into  a  low  pressure  part  and  a  high  pressure  part. 

line  of  the  operation  is  as  follows:  The  aqueous  feed  solution, 
•  s^a  water,  is  cooled  by  heat  exchange  with  out-going  fluid 
streams,  viz.  product  fresh  water  and  reject  concentrated  brine. 
In  order  to  balance  all  the  inefficiencies  of  the  process  plus  the 
srsible  work  for  separation,  a  refrigeration  step  would,  by  no- 
have  to  remove  a  certain  amount  of  heat  from  the  process. 
3  of  the  heat  should  be  removed  from  the  lowest  temperatu] 
-1  in  the  whole  process.   In  Fig.  ^,  it  is  shown  that  refri- 
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Fig.  3.   Flow    diagram    for   desalination    of     saline     water  by 
high    pressure    inversion     of    melting    points. 
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ration  is  used  to  remove  heat  from  the  recycle  organic  si  t 
id  bhe  aqueous  feed  solution  before  they  go  into  the  low  pres- 
sure converter  or  ice-maker.  T'ne   mixture  discharged  from  the  ice- 

:er  is  separated  into  an  organic  liquid  and  an  ice-brine  slurry 
in  separator  S^.   Ice-brine  slurry  is  separated  into  ice  and  con- 
centrated brine,  and  the  resulting  ice  is  washed  free  of  mother 
liquor  in  an  ice  washer.   The  ice  washing  operation  has  already 
been  well  developed  by  those  working  with  conventional  freezin 
operations  (2).   The  separated  concentrated  brine  is  heat  exch 
with  feed  solution  before  being  discharged.  Z!he   washed  ice  is 
reslurried  with  organic  liquid  to  form  an  ice  organic  liquid  slur 

1  the  mixture  is  pressurized  either  by  a  conventional  high  pres- 
sure pump,  or  by  a  flow-work  exchanger  (13)  •   '•••'hen  the  slurry  of 
ice  in  the  organic  liquid  is  pressurized  to  a  sufficiently  high 
pressure,  ice  begins  to  melt  and  the  organic  liquid  be  Ins  to  so- 

Fy. .  The  mixture  is  then  discharged  from  the  high  pressure  ice 
melter  and  is  separated  into  organic  slurry  and  fresh     r, 
are  e3    v   depressurized  through  conventional  turbines  or  by  flow- 
work  exchangers.   The  depressurized  slurry  is  further  coole 

ration,  as  has  been  previously  described,  and  is  ^     ;o 
.  ice-maker.   The  depressurized  fresh  water  is  ;.  _d 

.  ie  aq  „eous  feed  before  being  discharged  as  product 

_  *  thli       :  s .   I : .  iddit ion  ZQ   t . . 

to  all  freezing  processes,  th   .  .version  Pr       Proce: 

f Ol]  51 

1.   An  efficient  low-cost  orystalll 

due  to  the  condensed  yhuz^   operatic  . . 
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2.  Ice  crystal  formation  can  be  easily  controlled  because 
the  subcooling  can  be  kept  low  due  to  the  involvement 
of  only  condensed  phases.   This  advantage  leads  to  the 
possibility  of  large  granular  ice  crystal  formation  and 
thus  a  reduction  of  ice  crystal  washing  cost. 

3.  Energy  requirement  can  be  reduced  because  of  the  con- 
densed phase  operation  and  the  adoption  of  flow-work 
exchangers . 


IS 


EXPERIMENTAL 

One  of  the  objectives  of  this  work  was  to  set  up  a  continuous 

experimental  system  for  ice-making  at  atmospheric  pressure.   The 
ice  crystals  are  formed  according  to  the  following  equation: 

Aqueous  Solution  (C^)  +  Working  Medium  (solid) 


•Ice  +  Aqueous  Solution  (C2)  +  Working  Medium 

(liquid) 


(3) 


A  mixture  of  n-tridecane  and  n-tetradecane  was  chosen  as  the 

rking  medium  in  this  work  because  it  has  the  following  qualities: 

non-toxicity ,  a  very  low  solubility  in  water,  a  proper  melting 

range,  a  large  value  of  the  latent  heat  of  fusion,  a  high  value 

of  ( — )  for  melting,  no  hydrate  formation,  low  cost,  and  ready 
dp 

liability. 

;  ?ION  OF  IC   IAKER. 

-"ee  different  types  of  ice  maker  were  conside  ;ed, 

iey  were  1  (1)  a  counter-current  continuous-contact  vertical  eol- 
ith slow  agitation,  (2)  a  co-current  continuous-cont     to- 
rizontal  column  with  slow  agitation,  and  (3)  a  back-mi 
cter. 

[-mixed  type  contacter  was  finally 
>rk  for  the  following  reasonsi 
1.   '    nded  solids  can  be  h  idled  easily. 
The  d      of  dispersion  may  be 
thi         ileal  - 
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3.  Any  proportion  of  the  inlet  streams  aqueous  and  organic 
may  be  handled  easily,  with  either  phase  dispersed,  and 
the  ratio  of  the  two  flow  streams  may  be  controlled  at 
will. 

4.  Reasonably  reliable  scale-up  from  small  to  industrial 
size  is  possible, 

•"  II   "T  SET- UP. 

The  overall  experimental  set-up  which  was  used  for  obtaining 
the  data  is  shown  schematically  in  Fig.  5«   The  principal  piece 
of  equipment  used  was  the  ice-maker.   Other  items  of  the  equip- 
ment and  instruments  served  to  make  organic  slurry,  to  keep  the 
salt  water  at  its  freezing  temperature,  to  move  the  fluid,  to 
measure  the  flow  rate  and  to  achieve  the  recycling.   Figure  6 
shews  the  overall  view  of  the  equipment  used  in  the  closed  cycle 
operation.   A  description  of  the  individual  equipment  pieces  as- 
sociated with  the  fluid  flow  is  presented  below: 

Organic  slurry  maker.   The  equipment  was  modified  from  a 
soft  ice-cream  machine  made  by  the  Sweden  Electrical  Company. 
Figure  7  shows  the  schematic  diagram  of  this  slurry-making  system. 
Tr.e   exit  slurry  temperature  (and,  therefore,  the  solid  content) 
was  controlled  by  a  precision  temperature  controller  made  by  the 
Yellow  Spring  Instrumment  Company,  Inc.   It  was  an  "on-off"  type, 
si  71,  direct  dialing.   It  had  a  temperature  range  of  control 
from  -6°  C  to  -124°  C  with  a  sensitivity  better  than  ±0.05°  C 

slow  60°  C.   The  temperature  probe  was  installed  at  the  center 
of  the  cooling  cylinder. 


To  vacuum 


Permeability 


Manometer 


iO    vcc-~rr. 


V 
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Fig.    5.  Schematic   flow   diagram  of  a   closed    cycle  ice-making 
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Fig.  6.   Over  all  view  of  bench  scale  apparatus. 
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The  volume  of  the  cooling  cylinder  was  k   liters.   The  cooling 
rate  of  this  machine  was  about  200  cal/sec. 

Sslt  e  constant  temperature  bath .   This  bath  was  used  tc 

supply  the  salt  water  to  the  Ice  Maker  at  a  temperature  close  to 
its  freezing  temperature.   It  was  made  by  the  Blue  M  Electrical 
Co.,  Model  No.  PCC-A.   Temperature  was  maintained  by- a  cooling 
coil.   The  cooling  elements  were  controlled  by  a  "Blue  M  Microcon- 
troller".  The  capacity  of  the  bath  was  about  1^  liters. 

-nic  slurry  pump «   A  gear  pump  driven  by  a  C-.E.  1?25  rpm, 
one  horse  power  motor  was  used.   The  pump  was  made  by  Eastern 
Industries,  Model  No.  104-A21.  The   speed  of  the  pump  was  adjusted 
by  using  a  Zero-Max  torque  convertor,  Model  JKI,  speed  range  0  to 
^OOrpm.   The  maximum  capacity  of  this  pump  was  about  50  ml/sec. 

3 -.It  wat er  pump »      A  brass  gear  pump,  made  ^oy   the  Price  Pump 
Co.,  Model  No.  2,  driven  by  a  one  horse  power  motor  (Westinghouse, 
style  309?9l6A,  1?25  rpm)  was  used.   The  speed  of  this  pump  was 
regulated  by  using  a  Zero-Max  torque  convertor,  Model  1^2X,  speed 

nge  0-~00  rpm.   The  maximum  capacity  of  this  pump  was  about 
30  ml/  sec. 

Salt  water  flow  meter.   The  volumetric  flow  rate  of  the  salt 
er  from  the  Salt  Water  Pump  was  measured  cy   a  "Brook"  rota- 
:er,  Type  5-1110-5.   The  rotameter  was  calibrated  by  measuring 
-  time  required  to  fill  a  250  ml.  graduate  cylinder.   Calibra- 
>ns  for  the  temperature  of  3  C  were  run  and  are  given  in  Fig. 
B-l  in  Appendix  E. 

-  1c  slurry  temperature  recorder.   A  "  Honeywell  Two-Pen 
Electronic  19   Lab  Recorder"  was  used  to  indicate  and  to  record 


the  organic  slurry  inlet  temperature  of  the  Ice  Maker.   It  pro- 
vided high  speed  responses,  quick  selection  of  mul      s   n  (frc 
IOC  microvolts  to  100  volts),  positive  and  negative  zero  ran 
ability,  and  selection  of  multiple  chart  speeds  (10  fixed  speeds). 

pan  of   0.5  millivolts  was  used  throughout  the  experiment.   The 
accuracy  was  within  ±0.03°  C. 

~     ■'-■-.:?.      (See  Pig.  S)   The  salt  water  from  z'r.e   rotameter 
the  organic  slurry  from  the  pump  entered  the  Ice  Maker  from 
the  bottom.   The  vessel  was  constructed  of  §-in.  thick  plexl-glas . 
It  was  rectangular,  flat-bottomed  with  inside  dimensions  8.7  X 
3.?  X  16.5  cm.   The  agitator  impellers  were  attached  to  a  shaft 
located  at  the  center  of  the  vessel.   The  drive  motor  was  a  1/20 
HP  D.  C.  motor,  maximum  speed  1500  rpm.   The  agitator  speed  was 

julated  by  an  electronic  speed  controller  supplied  by  Cole- 
Parmer  Instrument  And  Equipment  Co.,  Catalog  No.  7025.  ^-e   re 
sed  was  measured  by  a  calibrated  tachometer,  obtained  from  S 
art-     c   Corporation,  Model  757-W.  To  minimize  the  heat  losses, 
the  outside  cf  the  vessel  was  insulated  with  1"  thick  styrofos 

From  the  Ice  Maker,  the  stream  was  directed  into  the  Recycle 
Vessel  to  achieve  the  recycling  or  into  the  Ice-  Sampling  Cc I 
to  obtain  samples  for  the  study  of  ice-crystal  size. 

V •   (See  Fig.  9)   This  vessel,  consist 

section  and  a  settler,  was  designed  to  make  po- 
of the  effluent  from  the  Ic       .  Is 
nic  solids  contained  In     effluent  frc 
ilted  in  the  heating  section,  then 
...  5  ..  sea  wore  separated  In 
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Fig.     8.    The    side  view    of    the    ice   maker. 
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2  was  returned  to  the  beginning  of  the  process. 
The  heating  section  was  agitated,  and  heat  was  supplied  by 
two  750-walt  stainless  steel  electric  immersion  heaters.   A  vari- 
able transformer  (made  by  The  Superior  Electrical  Co.,  Type 
2P71226,  0-2^-0  VAC,  2 A   KVA)  which  was  supplied  with  105  VAC,  was 
used  to  supply  the  desired  voltage  between  the  two  ends  of  the 
heater  A.   The  electric  power  to  the  heater  3  was  controlled  and 
regulated  by  an  "on-off"  type  temperature  controller  (made  by  the 
Yellow  Spring  Instrument  Co.,  Model  ?1.  For  detail,  see  the  tempe- 
rature controller  of  the  Organic  Slurry  Maker,  p.  19)  in  series 
with  a  variable  transformer  (made  by  The  Superior  Electric  Co., 
Type  216,  0.270  VAC,  0.81  KVA).   A  temperature  probe  was  installed 
at  the  exit  of  this  section  and  connected  to  the  controller.   By 
adjusting  the  output  voltage  of  the  transformers,  the  exit  stream 
temperature  could  be  kept  at  a  desired  value.   To  assure  complete 
melting  of  solids,  the  temperature  was  generally  controlled  at  +1 
C.   The  mechanical  agitation  in  this  section  was  accomplished  'oy 
two  marine-type  propellers.   The  drive  motor  was  supplied  oy   the 
Mixing  Equipment  Co. ,  Model  ?. 

The   phase  separation  in  the  settler  was  caused  by  the  differ- 

nce  in  the  density  of  the  two  phases  (0.78  gr/cm3  for  organic  and 

3  o 

1.02?  gr/cnr  for  aqueous  at  0  C).   The  volume  of  this  section  was 

rat  1.8  liters. 

""••'•-•anic  liquid  flow  i  eter.   The  flow  rate  of  the  organic 

stj     .  .z   measured  with  a  rotameter  after  it  left  the  settler. 

-  .  -scher-Porter  Flow  rater  tube,  No.  2F— :j- 25-5/70  was  used.  The 

rotameter  was  calibrated  by  measuring  the  time  required  to  fill  a 


25C  ml  graduate  cylinder  with  organic  liquid  at  0°  C.   The  j 
of  calibrations  are  given  in  Fig.  B-2  in  Appendix  5. 

___  : 1  Samrjlin.r  Column.   This  column  was  -     to  ocl"_ 

the  ice  crystals  in  the  effluent  of  the  Ice  Maker.   The  s     ty 
the  crystal  size  of  the  ice  was  then  determined.   This  colu] 
was  an  extruded  acrylic  pipe,  with  a  2- in.  ID,  2£-in.  OD,     1  - 
in.  long.   A  piece  of  60  mesh  screen  was  installed  at  t  -    attom 
of  the  column,  which  let  only  the  liquid  part  of  the  stres  -   iss 
through.  The  liquid  was  accumulated  in  a  closed  5  gallon  tank  to 

Lch  an  aspirator  was  connected  to  provide  the  drivi    "orce  of 
liquid  flow  from  the  column. 

A  permeability  measuring  column  was  installed,  but  was  n 
used  in  this  research  because  of  the  very  low  ice  production  ra 
of  the  system.   If  a  higher  capacity  system  is  studied,  the  perme- 
ability measuring  column  can  then  be  used  to  measure  the  sizes  cS 
ice   crystals  '  '-■'■■.  . 

-      _d  Insulation.   Due  to  the  strong  corrosiveness 

oho  salt  water,  Tyz^   tubing  was  used  mostly  in  this  syst<  .   rwo- 
incl    Lck  fiber  glass  was  used  to  insulate  the  outsi 
Qg  to  minimize  the  heat  losses. 

__;  -e  Salt  W ^  Throughout  r 

-  .  .01  Solution  was  used  instead  d~r   . 
[2]    that  the  pre 
atlon  very  similar  to  those  using  1 

.  -  ,d  by  disolvlng  table  salt  i. 


00 


ing  off  undisolved  material.  The  concentration  was  determined  by 

electrical  conductivity  measurements. 

?rcrr.ratlon  j   the  Or -rani  fixtures.   Mixtures  of  n-tridecane 
(d^Hos)  and  n-tetradecane  (C^^H^q)  were  used  as  the  working  med:  . 
in  this  work.   Their  physical  properties  are  as  follows: 

C13H28 
Molecular  wt.  1^8.25 

Melting  point  (1  atm. ) ,  (°C)  -6.2 

Eeat  of  fusion,  (cal/gram)   37  «0 

Solid  heat  capacity,         0.715(-10°C] 

(   r.o.1  /°n      rr-p   \ 


C14K30 

ref. 

198 

(15) 

-1-5.5 

(9  ) 

5^.4 

(15) 

0.418(-3°C) 

(15) 

Liquid  heat  capacity, 

(cal/°C  gr) 

Sp.  gr.  (20°CA°C) 


0.64l(-3°C)   0.515(-r7°c)   (15) 


0.757 


0.765 


(  9  ) 


Technical  grade  n-tridecane  and  n-tetradecane  were  obtained 
from  the  South  Hampton  Company,  Houston,  Texas.   Four  different 
kinds  of  mixtures  were  prepared  and  used  throughout  this  work, 
rheir  compositions  and  freezing  points  are  presented  belowi 


Mixture 
A 

■Q 

c 

D 


Composition    (Vol.) 
65%  n-C13H28,    35^  n-Cli;H30 
58>t  n-C13H28,    k2%  n-C-L^HoQ 
55%  n-C13H28.   ^5%  n-C^Hoo 
52%  n-C13H28,    L±&%  n-CliLj/H30 


Freezing  point,  WC 
-4>.^0  (measure 


-3.60 
-3.2c 
-3.07 


(measured) 
(measured) 
(measured) 


"ion  of  the  rocedure  for  an  Experimental  7/ 


Z -termination  of  Keat  Transfer  Hate 


(1)  The  Salt  Water  Constant  Temperature 
at   1°  G 

(2)  To  decrease  the  tine  required  for  the  sy     :o 
achieve  steady  state  temperature,  the  Salt  Water 
Pump  was  started  and  the  salt  water  (at  -1  C)  :■:. 
circulated  through  the  Ice  Maker  and  the  Recycle 
Column, 

(3)  The  Organic  Slurry  Maker  was  started.  TY.e   slurry 
temperature  was  controlled  at  a  desired  value. 

(s-)   The  Salt  Water  Pump  was  shut  off  temporarily  and 
the  Organic  Slurry  Pump  was  started.  The   torq 
convertor  was  adjusted  for  a  desired  organic  flow 
rate  as  determined  by  the  rotameter  a~  the  exit  of 
the  Settler.   Xeanwhere,  the  salt  water  exit  valve 
of  the  settler  was  closed. 

(5)  The  Salt  Water  Pump  was  started  again.  The  flc 
rate  was  adjusted  at  a  desired  value  as  det< 

oy   the  rotameter. 

(6)  The  agitator  of  the  Ice  Maker  was  started. 
rotary  speed  was  adjusted  as  determined  by  tl 
tachometer. 

(7)  The   agitator  of  the  Heating  Sect:. 

(8)  .  salt  water  valve  at  the  exit  of  the 
..  .  ■  julated  to  keep  the  J 

.idle  of  the  settl<  . 

(9)  The  power  to  t:       .-si on  he 

fter  the  tem]  of  tl 
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ing  Section  was  as  low  as  +1  C.   The  outlet  voltage 
of  the  transformers  was  adjusted  to  maintain  the 
temperature  at  -rl  C.   About  30  minutes  were  required 
for  conditions  to  reach  the  steady  state. 

(10)  The  experimental  conditions  were  recorded  on  a  data 
sheet  similar  to  the  one  given  as  Table  1. 

(11)  The  samples  of  brine  at  the  exit  of  the  Ice  Maker 
were  taken  through  the  sample  plug  by  using  a 
hypodermic  syringe. 

(12)  The  inlet  flows  to  the  Ice  Maker  and  the  agitator 
were  shut  off  simultaneously.   The  water  hold-up  in 
the  vessel  was  determined  after  the  ice  was  melted. 

3.   Determination  of  Salinity  and  the  Ice  Crystal  Size. 

(I)  (10)   The  same  as  those  in  Part  A. 

(II)  The  valves  to  the  Ice  Sampling  Column  and  the  He- 
cycle  Column  were  opened  and  closed  respectively  at 
the  same  time. 

(12)  Ice  crystals  were  accumulated  in  the  column.   The 
liquid  stream  was  directed  into  the  storage  tanki 
.  "-ermination  j   the  3rine  C one ent rat i en ,   An  electrical  con- 
ductivity meter  and  probe  were  found  to  be  a  quick  method  for  de- 

rmining  NaCl  content  of  the  sample  solution.   A  Model  RA4  Direct 
B  ading  Conductivity  Meter,  with  conductivity  range  between  0  and 
'■   microirihos/cm  and  with  a  probe  with  two  platinized  metal 
>s,  supplied  by  Industrial  Instruments,  Inc.,  were  used  for 
this  purpose.   This  concentration  measuring  system  was  calibrated 
reading  the  conductivity  of  NaCl  standard  solutions.      .-:- 


' 


Table   I.  .      sheet. 

Data  Shoot 

Run  No i  35  Date:  •#</?.  2&,   6  7 

__   slun 

coi  :ion ;    f&~ %  r- _C:13  flow  rate :  f.?5~  c .c./sec 

Ice-Maker  "^  enthalpy 

inlet  tamp.  \-4-S2     °C  of  fusion:  /^?    cal/sec 

Salt  water: 

composition:   3.5~       %  NaCl(wt.)     flow  rate:  3.5-    c. c./sec 
Ice-Maker 

inlet  temp.  :  -/*  j.  ^>       °C 

Agitation: 

speed:       3 o o  rpm  or      mech.   HPs  

Heater  power  input: 

Iff      amp.  X/gj"voit     or         ££/ cal/sec 

ing  column  exit  stream  temp.:  +0.5 °C 

3  — r  hold-up  in  Ice-Maker:      Soo         c.c. 

Ldence  time  of  salt  water  in  Ice-Maker:  /  43 sec 

.ner  liquor  analysis : 

sample:  1  2  3 

amt.   c.c.  /£  tj.         

vol.   of  dist. 

water  added:     2oo  £.00         

.iviLy 

>ter  ra    1        10.  <p_5~       /o.  9        

concentration 

( wt .        iCl ) :     44.29       4.3V       

sfer  from  salt  water: 

enthalpy  for  ice      s$.o 
nsible  heat:  Ho  cal/sec  crystal:-  1:  ss*t  cal 

total:    72,0.   7<sr. £    cal/sec 

efficiency:     jr/.g  ,    S3.  7 

tt: 
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suits  of  calibration  are  given  is  Pig.  3-3  in  appendix  3. 

To  analyze  a  salt  solution,  a  certain  amount  of  distilled 
water  was  added  to  dilute  it  so  that  its  concentration  -;:z.z   within 
the  meter  range  (600  to  1^00  ppm  NaCl).   The  meter  was  checked 
frequently  by  reading  the  standard  solutions  to  detect  any  change 
in  the  probe  or  meter. 

Ltat    ?oTA'er  Measurements.   The  cower  incut  to  the  im- 

—  -    ■■ — »■  ■■!■■!  HI   I-—.  ■  *  •*■ 

lers  during  the  operation  is  essential  for  the  scale-up  cf  the 
stem  (16).   In  this  work,  the  power  measurements  were  done   by 
measuring  the  electrical  power  input  to  the  drive  motor  with  an 
a;:meter  and  a  volt  meter  and  changed  to  horse  power  with  the  con- 
sion  1  h-p   =  7^5.5  watts,   A  correction  for  the  electrical  and 
chanical  power  losses  was  made  by  referring  to  the  motor  effi- 

Outout  cower  (agitation  ccwer) 


n^   out cut  cower  (agitation  r 
\       input  power  (measured  HP J 


Since  the  motor  manufacturer  was  unable  to  supply  the  infor- 
mation about  motor  efficiency,  data  for  various  load  conditions 
re  obtained  from  reference  books  (1?,  18). 

all  experiments,  air  was  removed  from  the  Ice-Maker  to 
prevent  the  aeration  of  the  liquid.  . 


RESULTS  AND  DISCUSSION 

^re  are  several  operating  variables  which  affect  the  ra 
of  productio.:  and  the  quality  of  ice  crystals.  Among  then  the 
following  variables  were  studied:  the  thermal  driving  force,  the 

;ree  of  agitation,  and  the  nominal  residence  time.  Experiments 
were  performed  to  investigate  each  of  these  while  the  other  vari- 
ables were  held  constant. 

The  thermal  driving  force  is  a  quantity  which  is  very  diffi- 
cult to  evaluate  in  this  complex  process.   Both  organic  and  aque- 
ous phases  change  their  composition  and  temperatures  simultaneous- 
ly and  therefore,  the  temperature  difference  between  phases  also 
changes  continuously.   In  fact,  the  driving  force  for  any  single 
small  portion  of  organic  and  aqueous  phases  can  be  quite  different 
from  any  other  meaning  that  the  thermal  driving  force  is  really  a 

itlstical  quantity.   However,  by  using  different  o:         tures 
with  different  incoming  temperatures  but  the  same  solid  content 

effect  of  the  thermal  driving  force  on  the  heat  transfer  rate 
could  be  studied  assuming  that  all  the  other  variables  were  ccn- 
.it.   5y  changing  the  salt  water  inlet  flow  rates  while  hold  . 
.  *t  organic  slurry  and  other  variables  constant,  the  effect 
of  thermal  driving  force  could  also  be  detected.   Ll'. 

he  temperature  and  composition  of  the  incoming  str 
..  and  varying  the  stirrer  speed  gave  the  effect  c 
t  transfer  rate.   Changing  the  total  inlet  fl 
the  !     flow  rate  ratio  between  two  incomin 
one  to  study  the  effect  cf  residence  time  or.  . 
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efficiency. 

The   recorded  experimental  data  are  summarized  and  shown  in 

Die  C-l  in  Appendix  C.   Calculations  for  a  run  were  made  as 
follows. 

First,  the  ice  crystal  content  of  the  product  stream  was 
calculated.   This  was  done  by  using  the  feed  and  product  assays 
of  the  salt  water.   Note  that  the  reaction. 

Aqueous  Solution  (C^)  -j-  Organic  Solid  —*»  ice  + 

Aqueous  Solution  (C2)  +  Organic  Liquid         [b) 

removes  fresh  water  from  the  salt  solution,  and,  therefore,  the 
salt  concentration  in  the  reactor  brine  increases.   Using  a  cer- 
tain amount  of  aqueous  feed  as  a  basis,  the  fraction  of  feed  water 
solidified,  and  hence  the  weight  percent s  of  ice,  salt,  and  water 
in  the  slurry  could  be  calculated.   Once  the  ice  content  was  de- 
termined, the  heat  transferred  for  the  ice  formation  could  be 
calculated  from  the  following  equation 

Qf  =  Wa  •  -rfo-  •  *i  (5) 

where  t 

Qf»  =  heat  transferred  for  ice  formation,  cal/sec. 

Wa  =  mass  flow  rate  of  inlet  salt  water  stream,  gr/sec. 

S  =  weight  percent  ice  in  product  stream,  %, 
Xj_  =  latent  heat  of  fusion  of  ice,  80  cal/gr. 

Next,  ~he  heat  transferred  for  the  cooling  of  salt  water 
from  the  inlet  temperature  to  the  freezing  temperature  was  calcu- 
I  in  a  straight  forward  manner  from  the  following  heat  balance 


Lon 

re: 

Qs  =  heat  transferred  for  cooling  of  salt  /rater  (or, 
sensible  heat  above  the  initial  freezing  point), 
cal/sec. 
t  n  a  =  salt  water  inlet  temperature,  °C. 


t*   ^  =  salt  water  initial  freezing  temperature,  -2.05  °C. 
C ^  a  =  mean  specific  heat  of  salt  water,  0.952  cal/gr,  °C 
for  3«5^  (wt«)  NaCl  solution  at  about  0°  C . 

The  total  heat  required  to  accomplish  a  known  decree  of  ic< 
'ormation  may  be  written  as 


+  Qf  +  Qu  +  Qi  (7) 


and  Qf  were   previously  defined.   Qu  and  Q^_  are  the  sensible 
it  of  unfrozen  brine  and  sensible  heat  of  ice,  respectively. 
Sin    he  temperature  of  the  product  stream  is  generally  not  lc 
-'.  -3*5  °C,  a  simple  calculation  may  show  that  the  sum  of  ~u 
account  for  only  a  small  portion  of  total  refr 
b.   Hence,  Qu  and  Q^  were  neglected,  and.  cense 

a  =  qs  +  ^ 

total  heat  transfer  efficiency  in  the  Ioe-] 
'olio-..  Li  -on 

6  =  — ! .  100 

Ao  *  ::f,o 


37 


re : 

3  =  total  heat  transfer  efficiency,  %. 
.. 0   =  organic  slurry  mass  flow  rate,  gr/sec. 
hf  0  as  heat  of  fusion  of  organic  solid  in  unit  mass  of  organic 
slurry,  cal/gr. 

The  calorinetric  determination  of  hf   0  and  calculations  of 
various  results  for  a  typical  run  are  presented  in  the  Appendix  C. 

Figure  10  is  a  plot  of  heat  transfer  rate  in  the  Ice-Maker 
versus  stirrer  speed  and  power  input  for  runs  No. 25  through  32. 
The  stirrer  speed  was  increased  in  increments  from  300  to  900  rpm. 
The  inlet  conditions  of  salt  water  and  organic  slurry  were  held 
constant  in  each  of  these  runs.   On  Figure  10,  the  topmost  curve 
is  the  total  amount  of  heat  available  for  transfer  from  the  heat 
of  fusion  of  organic  solid*   The  next  curve  down  shows  the  total 
amount  of  heat  transfer  per  unit  time,  which  was  used  in  the  for- 
mation of  ice  crystals  and  the  cooling  of  the  salt  water. 

lines  representing  the  enthalpy  -  temperature  relation  of  the 
jueous  system  (consisting  of  brine  and  ice)  and  that  of  the  organ- 
ic system  (consisting  of  solid  and  liquid  phases)  for  this  set  of 
n  is  are  shown  in  Fig.  11,  where  line  a-b-c  is  for  the  aqueous 
system  which  released  the  heat,  and  line  d-e-f  is  for  the  organic 

;cli  which  absorbed  the  heat.   Line  a-b-c  was  constructed  as 
follows.   Taking  the  salt  water  inlet  temperature,  t^n  a„  as  a 
reference  temperature  and  the  salt  water  mass  flow  rate  as  a  bas- 
.  the  enthalpy  of  the  aoueous  system  (H+.)  at  any  temperature  (t) 
jen  t^n  a  and  the  initial  freezing  temperature  of  salt  water 
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Fig.  10.     Heat    transfer   rate  vs.   stirrer   speed    and  po. 

sition  :  65%(vcl.)  n-C|3   ,  35%  n-C|4    ;    org.  slurry   f 

let  temp.  :-5.5°C  .    Salt  water  composition:  3-5%  (wt.)  ,    flow 

3.5  ml/sec   •  Inlet  temp.    ;  3  .G  °C  . 
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(tr   )  can  be  obtained  by  the  equation 
i  ,  a. 

Ht  =  Wa  -  Cp  •  (t  -  tln§a)  (1C) 

tfa  and  Cp  viere  as  previously  defined  for  Eq.  (6).   Line  a-b 
cc.    ponds  "Do  this  equation  which  represents  the  relation 

he  temperature  and  the  ethalpy  change  due  to  the  sensible  ..  sat 
evolution  of  the  aqueous  system.   At  temperature  tf   a,  ice  \ 
started  forming,  below  t  +   „,  (generally  designated  by  t*),  the  H- 
t  relation  can  be  written  as 

Ht,  =  [x  hlft,  +  (1  -  x)  hb§1..]  Wa  (11) 

.ere: 

.:  =  mass  fraction  of  ice  in  the  aqueous  system 
hj  i- 1  =  enthalpy  of  a  unit  mass  of  ice  at  t '  ,  cal/gr 
hv  +. »  =  enthaloy  of  a  unit  mass  of  unfrozen  brine  s~z   t •  , 
cal/gr. 

As  s 

hb,f  =  hlft.  -  Xi  (12) 

-  a_  is  the  latent  heat  of  fusion  of  ice,  Equation  (6 

"f  =  <hb,f  +  x*i)  '  -a  '   (13) 

z  •  c«      calculated  "::.         Jttion 

hbft.  =  cp  •   (f  -  ta  ti 

btained  by  usir. 
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(Fig.  A-4  In  Appendix  A)  and  the  mass  balance  relation.   A  rea- 
sonably linear  x-t  relation  was  obtained  for  the  operating  ran 

used  (as  shoi-m  in  Fig.  A-5 )  •   Line  b-c  was  constructed  by  usin 

is  figure  together  with  Equations  (13)  and  (14)  to  represent  the 
relation  between  temperature  and  the  enthalpy  change  after  ice 
crystals  started  forming. 

The  enthalpy-temperature  relations  of  the  organic  system, 
taking  the  slurry  inlet  temperature  as  a  reference  and  the  slurry 
mass  flow  rate  as  a  basis,  could  have  been  determined  in  the  same 
manner  as  that  of  the  aqueous  system.  However,  owing  to  the  lack 
of  a  suitable  phase" diagram,  some  critical  points  on  the  H-t  line 
were  determined  experimentally  instead.  On  Fig.  11,  point  "d" 
represents  the  condition  which  corresponds  to  the  inlet  slurry  tem« 

£ ;ure,  t*n  a,  and  zero  enthalpy  (tin  0  was  taken  as  reference 

temperature).   Point  "e"  was  found  hy   the  experimental  determina- 
tion of  the  initial  freezing  temperature  (t^  _)  of  the  specific 

J.  ,  u 

mixture  and  the  enthalpy  of  fusion  of  the  organic  solid  in  the 
slurry.   Because  of  the  lack  of  information  about  the  x-t  relation 
in  this  system,  which  was  required  for  the  construction  of  the  H-t 
relation  between  points  "d"  and  "e",  a  straight  line  was  drawn 
for  approximation.   As  the  temperature  went  higher  than  tf  0,  the 
enthalpy  change  with  temperature  was  completely  caused  by  'the  sen- 
sible  heat  absorption  and  was  determined  by  one  equation  similar 
to  Eq.  (10).   One  can  clearly  see  on  this  figure  that  the  salt 

stream  (at  3«5  ml/sec)  entered  the  Ice-Maker  at  +3.6°  C  and 
.s  cooled  to  about  -2  °C  after  19  calories  of  heat  had  been  trans- 
At  that  point  ice  crystals  started  forming.   One  also 


!■  o 


see  .  at  the  organic  slurry  entered  at  a  t€  srature  of  -5.5  °C 
an  lelted  absorbing  heat  from  the  salt  water.  The  format ic 
of  ice  continued  as  the  organic  solid  continued  to  melt.   For  a 

flow  rate  and  initial  temperature,  the  amount  of  contact  is 

nt  on  the  degree  of  mixing  in  the  contactor.        Id 

expect  an  increase  in  heat  transfer  rate  with  an  increasing 
of  mixing.  This  is  indeed  the  case  as  shown  by  the  ascending 
trend  of  the  curves  on  Fig.  10.   This  can  also  be  seen  by  comparing 
s  two  dashed  lines  parallel  to  the  vertical  axis  on  Fig.  11  show, 
ing  the  amount  of  heat  transferred  at  stirrer  speed  of  300  and  • 
.  respectively.-  It  should  be  noted  that  the  total  amount  of 
re  transfer  was  always  less  than  the  amount  of  heat  availabl . 
r  transfer  from  the  heat  of  fusion  of  the  organic  solid. 
.  sals  that  not  all  of  the  organic  solid  was  melted. 
Figures  12,  13 »  and  1^-  were  obtained  in  the  same 
Fig.  13  but  using  different  organic  mixtures  with  diff      Lnc 
:  -res  while  the  organic  solid  content  in  the  incomii 
slurry  in  ail  cases  was  about  the  same.   Taza  used  in  lonstructi: 
roe  figures  were  from  runs  No.  9-1?,  18-24,  and  ZZ-Z      -'  - 
ctive - 

From  t..    ta  used  in  FJ  i  .  10-14,  the  cu 
ncy  vs.  stirrer  speed  and  the  u>. 
F3  .15  and  Fig.  16,  respectively.   Due  :o  .  . 
tic  slurry  .       iure  and  th 
res  ill 
...    transf  • 
~sed  in  i 
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Fig.  12.     Heat    transfer    rate    vs.    stirrer      speed     and    power    input. 
Org.   mixture    composition  :   58  %  (vol.)   n  -  Cj2    7     42   % 
n-C,4     ;    org.    slurry    flow    rate  :  3.75   ml/sec    •,      inlet 
temp.    ;  -4.82°C.       Salt    water    inlet     composition:  3.5% 
(wt.)    NaCi    i     flow    rate     :    3.5  ml/sec    -,    temp.   •.  3.6°C. 
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Fig.  15.    HcQt  transfer   e  y    vs.  stirrei  er 
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.  .      force,   Here  again  ens  c  ..  note  that  the  total  ai 

ransfer  was  always  less  than  the  nt  of  con 

ich  could  be  accomplished  by  the  complete  melting  of  organic 
soli 

The  effect  of  the  thermal  driving  force  on  total  .. 
efficiency  can  also  be  evaluated  fron  another  set  of  sxperi- 
ntal  runs  (see  runs  No.  1-3  in  Appendix  c).   Mere  the  sal";  water 
inlet  flow  rate  was  changed.   However,  the  conditions  of  incoming 
organic  slurry  and  other  variables  were  kept  constant.   Figure  17 
illustrates  the  results.   Figure  18  shows  the  corresponding  tem- 
perature-heat  curves.  The   vertical  dashed  lines  indicate  the 
sunt  of  contact  in  each  of  these  runs.   One  can  see  that  the 
junt  of  heat  transfer  increased  even  though  the  mean  residence 
decreased.  The  decrease  of  the  mean  residence  time  was  causec 
by  the  increasing  salt  water  flow  rate.  The   increase  in  heat 
transfer  rate  with  the  increase  in  salt  water  flow  rate  can  best 

.  3x3 lained  by  the  increasing  salt  water  temperature  i    .e  pre  - 
uct  s      and  consequent  increase  in  the  driving  force  for  he 
jfer  as  shown  in  Fig.  18. 

linal  residence  time  is  defined  as  the  ratio  of  t 
vol       aqueous  slurry  in  the  Ice-Maker  to  the  vol 

salt  water.   The  effect  of  residence  time  on  th 

r  efficiency  is     ..  in  Pig.  19.  The  data  for 
c  tained  in  r    No.  ^0  thro   .  <7. 
1  flow  rates  were      3d  out  the  relativ 
..  and  0]         .  olid  co. . 
stirrer  s.   \    ..      pt  cor. 
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organic    mixture   composition:  55%  (vol.)   n-C 
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n-Ci4    ;   slurry   flow;  rate  :    8.75   ml     sec  •,   inlet   temperature 
4.52° c    .     Salt    wate  :sition  :  3.5  %  (wt.)  NaCI.    Stirrer 

id  :     400  rpm  . 
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of  residence  tine  on  the  amount  of  heat  transfer  could  also 

mined  by  ch   sing  contactor  size,   ho  attempt,  h    er,  was 
bo  do  so  because  it  would  require  a  change  in  the  si: 

so. 
rhe  ice  crystals  obtained  from  during  the  experiment 
studied  photographically.   Figures  20  and  21  she::  the  ice- 
crystal  pictures  which  were  taken  through  a  microscope.   As  can 
seen,  the  sizes  of  these  crystals  were  in  the  range  between  ~ . 3 
id  1.5  ram,  and  the  crystals  were  generally  rounded.   It  was 
found  that  thermal  driving  force  and  degree  of  agitation  did  if- 
fect  the  crystal  sizes.   A  smaller  thermal  driving  force  and  less 
agitation  resulted  in  larger  crystal  sizes,  and  vise  ver:  i.    ut 
.d  that  good  size  crystals  (larger  than  0.3 
0  could  be  obtained  under  almost  any  condition.   The  salini 
s  crystals  generally  ranged  bee     20  C  and  Z~ 
Lng  in  a  packed  bed  forced  in  the  Ice-Sam  li    '       ith  1; 
1.)  cold  fresh  water  to  the  amount  of  ice  obtained. 
call  th   .  Lgures  10  through  lo  she::  that  not  all 
organic  so]         Ited  in  each  run. 

•  used  was  of  the  back-n  I     .        herefor 
-  .  hi::  had  only  a  very  brief  c^  i 
L  left  3      ited. 

-         i      '  L me 

I< 

J  with 
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s  tracer.   The  salt  concentration  in  the  outlet  stream  was  meas- 
id  recorded  by  using  a  sensitive  conductivity  meter  and  a 
recorder,   roth  impulse  and  step  inputs  were  introduced  to  obtain 
the  C-curve  (a  response  to  an  upstream  impulse  signal)  and  F-curve 
(  i  response  to  an  upstream  step  signal). 

le  details  of  experimentation  and  data  analysis  of  this  part 
of  work  are  presented  in  Appendix  D. 

The   experimental  F-curves  and  C-curves  are  summarized  and 
ahc::r.   in  Figs.  22  and  23  which  are  the  same  as  Figs.  D-7   and  D-8 
respectively  in  Appendix  D.   The  corresponding  theoretical  curves 
for  various  extents  of  backmixing  -as  predicted  by  the  dispersion 
model  and  tank-in-series  model  are  included  for  comparison.   From 
these  two  figures,  one  can  see  that  the  behavior  of  the  Ice-Maker 
was  close  to  perfect  mixing  in  the  experimental  stirring  range. 

Using  the  experimental  C-curves,  the  parameters  which  charac- 
terize the  vessel  behavior  (i.e.,  extent  of  backmixing)  were  cal- 
culated.  The  results  also  indicate  an  approximately  "completely 
mixed"  condition  in  the  Ice  Maker. 

As  a  matter  of  fact,  the  residence  time  distribution  functions 
of  fluid  passing  through  the  Ice  Maker  as  found  from  tracer  studies 
are  only  partial  information  on  the  details  of  the  flow  patterns. 
that  is  known  is  the  length  of  time  certain  fractions  of  the 
fluid  stay  in  the  vessel  —  i.e.,  information  on  "macromixing"  (19) i 
n  :  the  molecule's  duration  of  stay  in  the  system  nor  Its 

th  other  molecules  which  it  encounters  during  its 
y   in  the  system  —  i.e.,  information  on  "micromi    "  (19*20) . 
In  general,  when  one  is  to  predict  the  performance  of  a  reac- 
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',  two  types  of  information  are  required!  (a)  the  f" 
rate  phenomena  (or  kinetics).   Since  age  distrj 
curves  produce  information  about  macromixing,  bat  not  micromixing, 
it  is  only  for  systems  in  which  the  interactions  between  iifj 

lecules  do  not  affect  the  rates  (i.e.,  first-order  rate  process- 
es) that  the  equipment  performance  can  be  predicted  solely  from 

~oer  curves  plus  kinectics  (21).  That  is  to  say,  the  performance 
of  a  reactor  is  independent  of  conditions  of  micromixing  (or  de- 
gree of  segregation)  if  the  rate  process  is  of  the  first  order. 
For  all  other  cases,  in  addition  to  age  distribution,  furt 
tails  about  the  mixing  are  required  in  predicting  the  performan 
of  equipment. 

In  this  work i  formation  of  ice  crystals  from  salt  water 
contacting  it  with  an  organic  slurry  was  ace. 

series  of  heat  and  mass  transfer  including  phase  cl        With 
her  complicated  process,  it  is  hard  to  conceive  Thai  t3 
-Z-  .  of  ice  is  a  first  order  process.   Theref     if  one 

o  sxplain  the  relationship  between  heat  t: 
.  speed,  as  shown  of  Figures  10  through  16, 
Lon  about  the  mixing  (i.e.,  micromi  I  -   or  sei 
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CONCLUSIC 

From  the  results  obtained  in  this  work  it  may  "be  concluded 
that  the  process  is  very  promising  because, 

(1)  The  size  of  the  ice-maker  will  be  relatively  small. 
From  the  rate  data,  it  appears  that  300  to  600  gallons 
of  fresh  water  can  be  obtained  from  each  cubic  foot  of 
the  contactor  per  day. 

(2)  The  ice  crystals  produced  in  this  process  were  fairly 
large,  well  shaped,  smooth  and  low  in  salt  content. 
This  indicates  that  the  ice-washing  operation  can  be 
easy  to  carry  out. 


PROPOSED  E     I  OF  A  SEMI-PILOT  PLANT 

ta  obtained  in  this  work  together  with  thermodynamic  lata 
the  working  medium  obtained  by  S.  3.  Maa  car;  be  u    1  "or 
designing  a  continuously  operating   1  jrfc. 

Figures  2~  and  25  show  the  tentative  design  for  a  thousai 
lion  per  day  pilot  plant.   Table  2  gives  the  nomenclature  and 
sizes  of  equipment  shown  in  Figures  2^  and  25. 

3  major  piece  of  equipment  that  has  not  yet  been  investiga- 
ted is  she  ice-washer.   But  the  design  and  operation  of  this 

at  have  already  been  well  developed  in  connection  with  the  con- 
ventional freezing  processes.  Information  could  be  obtained  from 
literature. 

Dther  major  piece  of  equipment  not  yet  st     1   is  the 

ice-melter.    -t  since  ice-melting  is  simpl;   lie  reverse 
operation  of  ice-:..    j,  it  can  be  expected  that  the  heat  transfe] 
rats  will  be  similar  to  that  in  the  ice-maker. 
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APPENDIX  A 

PHYSICAL  PROPERTIES  0? 
WORKING  MEDIUMS  AND  SALT  WATER 
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CALIBRATIONS  0?  INSTRUMENTS 
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Fig.  B-l.  Salt  water  flow  meter  calibration  , +3°C   . 
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rhe  experimental  data  are  listed,  in  Table  C-3»  The  necessary 
sulations  are  presented  below. 

1.   Calculation  of  heat  transfer  rate  for  ice  formation. 
From  sar.pl e  !•• 

4  ml  of  mother  liquor  was  diluted  with  200  ml  of 
distilled  water.   Conductivity  measurement  showed  the 
diluted  solution  has  a  reading  of  10.85  on  the  meter, 
which  corresponds  to  a  NaCl  cone,  of  842  ppm. 

Hence,  using  mass  balance,  wt .  %   of  NaCl  in  the 
mother  liquour  was 


204      842      ..     ,, 
IT  X  1,000,000  x  10°  =  4'29 


In  the  feed  stream,  the  wt.  %   of  NaCl  was  3»5- 
Taking  100  gr.  of  feed  as  basis,  and,  again,  using 
mass  balance,  the  mass  of  ice  formed  in  the  product 
stream,  S,  was  estimated  from  the  equation 

100  x  0.035  =  (100  -  S)  X  0.0429 

This  gave 

S  =  18.4  gr/100  z'r   of  feed. 

That  is,  in  100  gr  of  salt  water  feed,  18.4  gr  of  ice 
3  produced. 

The  salt  water  volumetric  inlet  flow  rate  was  3»5 
ml/sec  corresponding  to  the  mass  flow  rate  of 


Wa  =  3.5  x  1.02?  -  3.595  gr/sec. 
lerefore,  ice  production  rate  was 


3.595  x  t§~  =  °-661  gr/sec. 


Finally,  the  hear  transfer  rate  for  ice  formation 
wa s  cc mp ut  ed  as 

Qf  =  0.661  x  80  =  53.0  cal/sec. 

Prom  sample  2: 

Again,  4  ml  of  -other  liquor  was  dilute!  with  .. 
ml  of  distilled  water.   The  conductivity  measurement 

showed  the  diluted  solution  had  a  concentration  of  J 
trom  and  thus 


*  ~  4  X  1,000,000  x  ±uu  -  *•**• 


Further  more,  the  material  balance, 

100  x  0.035  =  (100  -  S)  x  0.0^-34, 


3  =  19.35  gr/100  sr  of  feed. 
;refore,  the  ice  production  rate  in  this  ru: 


3.595  x  — --  =  0.695  gr/:   , 


.  for  ice  f 
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Qf  =  0.695  x  80  =  55.6   cal/sec. 

2.   Calculation  of  the  heat  transfer  rate  for  the  cooling  of 
salt  water. 

t1r  a  =  3.6°  C   • 
in  >  a. 


tf  0  =  -2.05  C 


C-  H  =  0.952  cal/gr  °C 

Wa  =  3.595  gr/see 

From  equation  (6),  Qs  was  calculated  as 

Qs  =  3.595  [3.6  -  (-2.05)]  0.952 

=  19  cal/sec. 

3*   Calculation  of  total  heat  transfer  rate. 
From  sample  li 

Using  equation  (8), 

Q  =  Qf  +  Qs  =  53.0  +  19.0  =  72.0  cal/sec. 

Prom  sample  2s 

C;  =  Q^  +  Qs  =  55«&   +  19.0  =  7^.6  cal/sec. 

k.      Calculation  of  total  heat  transfer  efficiency. 

The  hear  of  fusion  of  organic  solid  entering  the  Ice- 

r  T.:as  139  cal/sec. 
Ac c o r d i ng  to  E quat ion  ( 9 ) 


Q   =  -.-.-■  =  .  ._  ~   or  5:  • 


3  =  __  =  0.537  or  53.7^        pie  2. 


5,   Calculation  of       m  reside 

■    -       lurry  >00  ml 

and  the  volumetric  inlet  flow  rate  of  s  ;  3.5 

/sec. 

is,    u>o  no 1  2  .  .  was  -^4  =  I 
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Determination  of  hf ,o,   The  heat  of  fusion  of  organic  solid 
was  determined  calorimetrically .   A  200  ml  vacuum  bottle  provided 
with  a  cork,  a  ring  stirrer  made  of  copper  wire,  and  an  electric 
immersion  heater  with  44.5  ohms  resistance  was  used  as  a  calorim- 
eter.  An  ammeter  with  an  accuracy  of  0.1  ampere  and  a  variable 
transformer  were  used  to  measure  and  adjust  the  current  to  the 
heater.   The  temperature  recorder  described  on  page  23  was  used 
to  indicate  and  record  the  temperature  in  the  bottle. 

The  total  energy  required  to  heat  the  whole  system  from  its 
initial  temperature  to  a  temperature  a  few  degrees  Centigrade 
higher  than  the  freezing  point  of  the  organic  mixture  could  be 
written  as 

K  =  Hs  +  Hf>Q  +  Hu  +  Hi  +  Hc  (o-l) 

where: 

H  =  total  energy  required  for  heating,  cal. 
Hs  =  sensible  heat  of  organic  liquid  above  the  initial 
freezing  temperature,  cal. 
Hf   =  heat  of  fusion  of  organic  solid,  cal. 

Hu  =  sensible  heat  of  organic  liquid  below  the  initial 

freezing  point,  cal. 
EL  =  sensible  heat  of  organic  solid,  cal. 
K.  =  heat  absorbed  by  the  calorimeter,  cal. 

Hs  and  Hc  could  be  calculated  from  the  following  equations 
Hs  =  Mo  *  Cte  -  tffQ)  •  Cp>ol  (o-2) 
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Ho  =  ^e  "  *lni>  *  CP.c  <°-3) 

where: 

!•:  =  mass  of  organic  mixture  in  the  calorimeter,  gr. 

t.  .  =  initial  slurry  temperature,  °C . 

mi 

t  =  final  temperature  of  heating,  °C . 
t~   =  organic  mixture  initial  freezing  temperature,  °C . 

1,0 

C^  0]_  =  specific  heat  of  organic  liquid,  cal  deg"1  gr"1. 
Cp  c  =  heat  capacity  of  calorimeter,  cal  deg"1. 

cp  02_  was  determined  by  the  composition  of  the  mixture  and 
specific  heat  of  n-tridecane  and  n-tetradecane.   C^  r   was  deter- 
mined  experimentally. 

H„  and  K^  in  Eq.  (c-1)  were  neglected  because  they  account 
for  only  a  very  small  portion  of  K. 

Hence,  Equation  (c-1)  becomes 

H  =  MQ  •  (te  -  tf >0)  •  Cpt0l  +  Hf >0  +  (te  -tlnl)  •  CpfC 

• (c-4). 

Z!he   heat  absorbed  by  the  system  equals  to  the  total  heat 
input  because  of  the  conservation  of  energy. 

Neglecting  the  energy  input  from  the  stirrer  and  surround- 
ings, the  total  energy  input  was  from  the  electric  heater,  that 
is, 

2 

E  =  i2  Rt  joules  =  ltlH   cal  (c-5) 

>re « 

E  B  total  electrical  energy  input,  cal. 
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i  =  electric  current  to  the  heater,  amp. 

R  =  heater  resistance,  ohms. 

t  =  total  heating  time,  seconds. 

Equating  Eqs.  (c-4)  and  (c-5)  yielded 

Hf,o  -  klM   *  M°(te  "  tf.o)  '  CP,ol  -  (*e  -  tini)  '  CPfl 


(c-6). 


The  heat  of  fusion  of  organic  solid  in  a  unit  mass  of 
organic  slurry  could  be  obtained  by  dividing  both  sides  of  Eq. 
(c-6)  with  M0. 

hf  n   =  -±*£ 
1  » o    M0 

i   Rt  1 

=  4.184  M0  "  (te  "  tf,o)  '  cp,ol  "  m;  (te  "  tini),Cpt 

The  results  of  this  part  of  work  are  given  below. 
Mixture     Composition  (vol.)     Initial  Slurry     hf  0>cal/gr 

Temperature,  tinj_   °C 


65fo   n-C13H2s 
35%   n-ClZ|H30 


-5.50  17.73 


55%  n-C13H28 
3  J  -4.82  29.73 


k2%  n-ClLiE30 

55%  n-C13H28 
^5%  n-ClLiZ30 


-4.52  20.37 
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52%  ^-Ci3H23 
48*  n-ClZfH30 


-4.35 


20.95 
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APPENDIX  D 

ANALYSIS  OF  RESIDENCE  TIME  DISTRIBUTION  OF  FLUID  IN 
ICE  MAKER 
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Two  ideal  flow  patterns  which  are  often  used  to  approximate 

real  systems  are  plug  flow  and  complete  mixing.  Patterns  of  flow 
other  than  plug  flow  or  complete  mixing  are  called  non- ideal  flow 
patterns  (22) . 

Knowledge  of  the  complete  flow  pattern  of  the  fliid  In  the 
system  is  required  if  one  is  to  account  for  non-ideal  flow.   Due 
to  the  difficulties  associated  with  obtaining  and  interpreting 
such  information,  an  alternate  approach  is  used  which  requires 
knowledge  only  of  how  long  different  elements  of  fluid  reside  in 
the  system.   This  partial  information  is  relatively  simple  to  ob- 
tain experimentally,  can  be  easily  interpreted,  and  yields  infor- 
mation which  in  many  cases  allows  a  satisfactory  accounting  of 
the  non-ideal  flow. 

The  experimental  technique  used  for  finding  this  desired 
distribution  of  residence  times  of  fluid  elements  in  the  system 
is  a  stimulus  -response  technique  using  some  sort  of  tracer  mate- 
rial In  the  inlet  fluid  stream.   The  injection  is  the  stimulus 
and  the  response  is  the  tracer  concentration  measured  in  the  out- 
let stream.   The  tracer  can  be  radio-active,  or  a  colored  dye,  or 
a  salt  solution,  or  any  material  depending  on  the  particular  situ- 
ation.  Tracer  is  injected  into  the  inlet  stream  in  some  known 
fashion,  such  as  a  step,  a  pulse,  a  semi-wave  or  other  cyclic  sig- 
nal, or  even  a  random  signal  with  known  properties. 

Quantitative  definitions  and  discussions  of  the  age  distr 
tion  functions  developed  by  Danckwerts  and  others  could  be  f 
In  the  literatures  (22,23).   Some  particular  terms  used  in  this 
work  are  discussed  below i 
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The  ?-Curve.   With  no  tracer  initially  present,  let  a  step 
tracer  signal  of  concentration  C0  be  introduced  into  the  fluid 
entering  the  system  in  such  a  manner  that  the  volumetric  flow  rate 
to  the  vessel  remains  constant.   Then  the  concentration-time  curve 
for  tracer  in  the  exit  fluid  stream,  measured  in  terms  of  inlet 
.tracer  concentration  and  reduced  time  is  called  the  F-Curve.  The 
range  of  F  is  0  to  1. 

Reduced  time,  6,  is  defined  as 


(D-l) 


where  t,  the  mean  residence  time  of  fluid  in  a  vessel  is  defined 
as 

t-  _    Volume  of  the  vessel  available  for  flow (T)-2) 

Volumetric  flow  rate  of  fluid  through  the  vessel   ~ 

Thj  C-Curve.   The  curve  which  describes  the  concentration- 
time  function  of  tracer  in  the  exit  stream  of  a  system  in  response 
to  a  delta  function  or  unit  impulse  input  is  called  the  C-Curve. 
As  with  the  F-Curve,  the  range  and  domain  are  dimensionless. 
Concentrations  are  mearsured  in  terms  of  the  initial  concentration, 
C°,  as  if  it  were  evenly  distributed  through  out  the  system. 

models  FOR  NONIDSAL  FLOW 

There  are  many  types  of  models  which  can  be  used  to  charac- 
terize nonideal  flow  patterns  within  vessels  (2*0.   Two  models 
which  will  be  employed  are  briefly  discussed  below i 

Dispersion  Model  (Dispersed  Plug  Flow  Model) .   This  model  was 
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drawn  on  the  basis  of  analogy  between  mixing  in  actual  flow  and  a 
diffusional  process,  with  a  constant  dispersion  coefficient  and  a 
flat  velocity  profile  in  the  vessel  considered.  The  parameter 
which  characterizes  the  degree  of  back  mixing  during  flow  is  a 
dimensionless  group  D//J   L,  called  the  vessel  or  reactor  dispersion 
number.   It  changes  from  zero  for  plug  flow  to  infinity  for  back- 
mix  flow  and  is  the  reciprocal  of  the  axial  Peclet  number  for  mass 
transfer. 

To  characterize  flow  in  a  real  vessel,  one  needs  to  select 
from  the  family  of  theoretical  F  or  C  curves  the  curve  that 
most  closely  fits  the  experimental  F  or  C  curve.   From  the  D/jjl 
of  this  fitted  curve  the  behavior  of  the  vessel  as  a  reactor  can 
then  be  predicted. 

Another  convenient  way  to  predict  the  reactor  behavior  is  to 
equate  the  variance  of  an  experimental  C  curve  to  that  of  a  se- 
lected theoretical  C  curve.   For  the  C  curve  in  a  closed  vessel, 

&         i   _L  .  o    <JL\2    (l-v-ML/D)  (d-3) 


<^§=^-2(*r<l-e-^/D> 


where 


t  =  mean  age  of  the  exit  stream 

2       ? 
£q,  and  tft  =  variance  of  the  C  curve,  which  measures  the 

spread  of  the  distribution  about  the  mean. 

From  Sq.  D-3»  nr  for  a  actual  reactor  could  be  solved  for 

2 

with  a  tfQ  obtained  from  the  experimental  C-curve. 

Completely  Stirred  Tank  '.       3  -srles  Model,   This  model  is 
based  on  the  assumption  that  the  actual  reactor  can  be  repr 
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by  a  series  of  j  equal  sized  backmix  flow' vessels.   The  output  C 
curve  for  this  model  is 


C  =  KSu       *'**  (D-M 


with  mean      9C  =  1 


and  variance   5Q  =  ■— 


From  experimental  variance  measurements  j  can  be  found. 


EXPERIMENTAL. 


In  this  work,  the  residence  time  distribution  of  fluid  in 
the  Ice  Maker  was  measured  with  a  steady  flow  of  fresh  water  and 
an  injection  of  saturated  sodium  chloride  solution  as  a  tracer. 
A  10  ml.  syringe  was  used  to  create  an  impulse  input,  and  a 
"Compact  Infusion  Pump",  made  by  Harvard  Apparatus  Company,  Model 
975.  with  a  50  ml.  syringe  to  create  a  step  input.   A  conductivity 
cell  was  installed  at  the  outlet  of  the  vessel  along  with  a  Model 
RA4  Direct  Reading  Conductivity  Meter  (see  p. 31)  a^d  a  Honeywell 
Electronic  19"""  Lab  Recorder  (see  p.  23)  which  gave  a  record  that 
was  both  continuous  and  linear  with  concentration. 

The  experimental  runs  were  made  with  stirrer  speeds  of  300, 
600,  and  900  rpm. 

RESULTS  AMD  ANALYSES. 

The  experimental  conditions  for  each  run  are  summarized  and 
shown  in  Table  D-l. 


Table  D-l 

Summary  of  experimental  conditions  for 
each  run  of  residence  time  distribution  study. 


Hun  Xo. 

Type  of 
input 

Water  flow 

rate 

ml /sec 

Mean  residence 
time,  sec.+ 

Stirrer 
speed 

1 

step* 

12.1 

103 

3C0 

2 

step"* 

12.1 

103 

600 

3 

step* 

12.1 

103 

900 

it 

impulse* 

12.1 

103 

300 

5 

impulse""- 

12.1 

103 

600 

/■ 

imDulse* 

12.1 

103 

900 

+  vessel  volume  12^+9  ml, 

*  the  tracer  injection  rate  was 
0.103  ml/sec. 

*  the  amount  of  tracer  used  was  10 
ml.,  time  width  about  3  seconds. 
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The  experimental  traces  of  individual  runs  are  reproduced  in 
Figs.  D-l  through  D-6.   The  recorded  curves  are  quite  smooth. 

The  tracer  concentration  readings  obtained  from  the  exper- 
imental traces  are  presented  in  Tables  D-2  and  D— 4.   Tables  D-3, 
D-5 »  a^d  D-6  show  the  calculated  values  of  F(6),  D/uL   and  j,  and 
C(9)  respectively.   Figures  D-l  and  D-2  are  the  plots  of  F-curves 
and  C-curves.   The  theoretical  curves  for  various  extents  of  back- 
mixing  as  predicted  by  the  dispersion  model  and  tank-in-series 
model  are  included  for  comparison. 

In  comparing  the  experimental  F-curve  with  theoretical  ones, 
one  can  see  that  the  behavior  of  the  Ice-Kaker  in  the  operating 
stirrer  speed  range  was  very  close  to  the  "  completely  mixed"  con- 
dition.  But,  the  C-curves  and  the  calculated  values  of  mixing 
parameters  indicate  that  there  was  a  small  degree  of  deviation 
from  the  conditions  as  F-curves  implied.   This  was  due  to  the 
assumption  of  "impulse  input"  in  the  analysis.   Actually,  the  true 
impulse  should  be  done  by  injecting  an  finite  amount  of  tracer 
over  an  infinitesimal  period  of  time,  and  this  is  really  a  physi- 
cal impossibility.   The  assumption  of  tracer  injection  in  this 
work  (injection  time  width  about  3  seconds,  mean  residence  time 
of  fluid  in  the  vessel  about  100  seconds)  is,  in  fact,  an  approx- 
imation.  Therefore,  it  may  be  concluded  that  the  F-curve  gives 
a  better  interpretation  of  reality,  and  the  behavior  of  the  vessel 
was  close  to  perfect  mixing   (or  well-stirred). 


IOC 


PHIWOOyHIA,  PA 

-r— 


PRIi.TED   IN  U.  S.  A. 


O  I     N^/_C  ' 


Fig.  D-l.   Response  curve  of  run  No.l  (step-input, 
300rpm) 
Recorder  chart  speed i  0.5  in/mi n. 
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Fig.  D-3.   Response  curve  of  run  No. 3  (step-input,  900  rpm) 
Recorder  chart  speed:  0.5  in/rain. 
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Table  D-2 
Experimental  data  (step-input  response) 


Sun 

No.  1 

Hun 

No.  Z 

Him 

WO.  j 

Time 
( sec 

:  Tracer  Cone. 
. )     (ppm) 

Time  Tracer  Cone, 
(sec)     (ppm) 

Time 
( sec 

:  Tracer  Cone. 
0    (ppm) 

0 

0 

0 

0 

0 

0 

20 

260 

30 

370 

24 

310 

50 

550 

60 

£50 

54 

620 

30 

765 

90 

865 

84 

830 

110 

940 

120 

1050 

114 

1010 

140 

1080 

150 

1175 

144 

1145 

170 

1180 

180 

1270 

174 

1280 

200 

1265 

210 

1340 

204 

1320 

230 

1320 

240 

1420 

234 

1330 

260 

1380 

270 

1455 

264 

1430 

290 

1450 

300 

1490 

294 

1455 

320 

1465 

330 

1520 

324 

1490 

350 

1480 

360 

1530 

354 

1510 

330 

1435 

390 

1545 

334 

153C 

410 

1500 

420 

1550 

414 

1540 

440 

1500 

450 

1550 

444 
474 

1550 
1550 

L/    ZZ 

103 

sec 

0  = 

103 

sec 

~  m 

103 

sec 

Co  B 

1500 

pm 

Co  = 

155C 

1  ppm 

Co  = 

1550 

D 

10? 


Table  D-3 
Dimensionless  data  (step-input  response) 


e 

t 

F(9)'  =  -£i- 

c0 

Run 

No.    1 

Run 

No.    2 

Run 

No.    3 

0 

F(0) 

e 

F(G) 

8 

7(9) 

0 

0 

0 

0 

0 

0 

0.194 

0.173 

0.291 

0.239 

0.233 

0.20 

0.485 

O.367 

0.583 

0.419 

0.524 

0.40 

0.777 

O.510 

0,.874 

0.558 

0.815 

0.535 

1.067 

O.626 

1.165 

0.677 

1.107 

0.652 

1.359 

0.720 

1.456 

0.758 

1.398 

0.?38 

1.650 

0.787 

1.7^8 

0.819 

1.689 

0.825 

1.942 

0.843 

•     2.039 

O.865 

1.981 

0.852 

2.233 

0.880 

2.330 

0.916 

2.272 

0.890 

Table  D-4 

Experimental  data  ( impulse- input  response) 
(smoothed  to  equidistant  points) 


Run  No. 4 

Run  No. 5 

Run  No. 6 

(300  rpm) 

(600  rpm) 

(900  rpm) 

Time  (sec) 

Concentration 
(ddih  NaCl) 

Concentration 
(■opm  NaCl) 

Concentration 
(r;r~.    NaCl) 

0 

0 

0 

0 

12 

740 

750 

700 

24 

810 

800 

743 

36 

740 

750 

693 

48 

675 

685 

630 

60 

615 

625 

560 

72 

555 

565 

510 

84 

500 

505 

460 

96 

450 

450 

420 

108 

400 

410 

370 

120 

360 

370 

330 

132 

323 

323 

300 

144 

235 

290 

260 

156 

250 

255 

235 

168 

222 

225 

205 

180 

200 

205 

180 

192 

173 

180 

160 

204 

152 

160 

140 

216 

135 

140 

122 

228 

120 

125 

105 

109 


Table  D-4( Cont'd) 

Experimental  data  (impulse-input  response) 
(smoothed  to  equidistant  points) 


Run  No.  4 

Run  No. 5 

Run  No. 6 

(300  rpm) 

(600  rpm) 

(900  rpm) 

Time  (sec) 

Concentration 

Concentration 

Concentration 

(ppm  NaCl) 

(ppm  NaCl) 

(pDm  NaCl) 

240 

105 

110 

95 

252 

90 

100 

85 

264 

80 

88 

75 

276 

70 

75 

65 

288 

60. 

65 

55 

300 

50 

55 

48 

312 

43 

50 

40 

324 

36 

45 

35 

336 

30 

35 

30 

348 

25 

25 

25 

360 

20 

20 

20 

110 


Table  D-5 
Calculation  Results  (impulse-input  response) 


Hun  No.  4 

Run  "o. 

5 

F  n  No,  6 

IC;L 

8314 

8481 

.7696 

ECjAt  =  12  x  2C1 

99,768 

101,722 

92,352 

SC^i 

824,196 

853,968 

762,372 

"   ZCi   (sec) 

99.13 

100.69 

99.06 

1.  1 

130,429,872 

137,313, 

,400 

■  121,147,488 

^C ,  t ,     _  ^ 

^  2            11         /  X.  \ 

5861.22 
0.5965 

6052.82 
0.5970 

6,413.68 
0.6536 

6t   -  2Cl  "  <*' 

0©  -   -x 

(t) 

Dispersion  Model 
Parameter  (D^t/L)       0.557         0.558         0.698 

Tank  in  Series  Kodel 
Parameter  (j)        I.6766        1.675         1.5299 


Ill 


Table  D-6 
Diraensionless  data  ( impulse- input  response) 

G  =  t/  t 

Run  No.  4             Run  No.  5             Run  No.  6 

9   G(  9  ) 9  Ci(9) 9  C,-(  e  ) 

0           0                                           0           0                                          0  0 

0.121           0.735                       0.119           0.742  0.121  0.751 

0.242          0.805                     0.238          0.791  0.242  0.797 

0.363          0.735                     0.358          0.742  O.363  0.743 

0.484         0.671                    0.477          0.678  0.484  0.676 

0.605          0.611                    O.596         0.618  0.605  0.601 

0.726          0.551                    0.715         0.559  0.726  0.547 

0.347          0.497                    0.834         0.500  0.848  O.493 

0.968          0.447                    0.954         0.445  0.969  0.450 

1.089          0.397                     1.073          0.406  1.090  0.397 

1.210          0.358                     1.192          O.366  1.211  0.354 

1.331        0.321                 1.311        0.320-  1.333  0.322 

1.453          0.283                     1.430          0.287  1.454  0.279 

1.574          0.248                     1.549          0.252  1.575  0.252 

1.695           0.221                      1.668           0.223  I.696  0.220 

1.816          0.199                    1.788          0.203  1.817  0.193 

1.93?          0.172                    1.907          0.178  1.938  0.172 

2.058          0.151                    2.026         0.158  2.059  0.150 
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Fig.  D -8.  Impulse    input    response   (or  C-)   curves. 
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ABSTRACT ' 

Experimental  equipment  was  designed  and  built  to  test  the 
ice-making  of  the  Inversion  Desalination  Freezing  Process  at  at- 
mospheric pressure.   A  back-mixed  type  contactor  was  used  to  pro- 
duce ice  from  salt  water  by  direct  contact  freezing  with  a  mixture 
of  n-tridecane  and  n-tetradecane  as  a  working  medium. 

The  thermal  driving  force,  degree  of  mixing,  and  nominal 
residence  time  were  considered  as  controlling  variables  for  the 
ice  production  rate  and  have  "oeen   studied.   Forty  seven  experimen- 
tal runs  were  carried  out  using  four  different  working  media. 
The  agitation  speed  ranged  from  300  to  900  rpm,  and  the  nominal 
residence  time  from  86  to  238  seconds.   In  addition  to  the  ice 
production  rate,  the  sizes  of  ice-crystals  obtained  under  differ- 
ent operating  conditions  were  investigated  photographically.   The 
pictures  showed  the  crystals  were  fairly  large,  well  shaped,  and 
smooth. 

The  results  of  the  production  rate  study  can  be  used  to  design 
a  larger  scale  plant.   The  experimental  data  reveal  that  very 
reasonably  sized  equipment  may  be  used  with  this  process.   The 
quality  of  ice  indicates  that  the  ice-washing  should  present  no 
'or  obi  ems. 


